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ABSTRACT 

Analysis of uncertainties in stratospheric perturbations 
requires the development of a model which uses a mini- 
mum of computer time yet is complete enough to reason- 
ably repi-esent the results of more complex models. 

Such a concise model has been developed and is described 
in this paper. The model is a steady-state model using 
iteration to achieve coupling between interacting species. 
The species NjO, odd nitrogen, CH^ , CH^ Cl, CCl^ , 
CFjClj, CFClj, and odd chlorine are determined from 
diffusion equations with appropriate sources and sinks. 
The HOj, species, 0, 0*D and are solved in photo- 
chemical equilibrium and H 2 C, I-I 2 , N 2 , O 2 and temper- 
ature have been assigned fixed profiled. Diurnal effects 
due to chlorine nitrate formation are accounted for by an 
analytic approximation. The model has been used to eval- 
uate steady-state perturbations due to injections of chlo- 
rine and NOjj. The results are similar to those obtained 
by other models. 
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MONTE-CARLO ANALYSIS OF UNCERTAINTY PROPAGATION 
IN A STRATOSPHERIC MODEL: 

I. DEVELOPMENT OF A CONCISE STRATOSPHERIC MODEL 


INTRODUCTION 


In the study of the potential effects of man-made pollutants on the Earth's 
ozone layer, computer models of the stratosphere have been used extensively 
to provide quantitative predictions of ozone depletion. These models utilize 
as inputs experimentally determined transport coefficients, chemical reaction 
rates, and photodissociation cross-sections, and produce as outputs the con- 
centrations of various chemical species as a fur-;tion of altitude. In order to 
make an intelligent assessment of the impact of stratospheric pollution, it is 
necessary to know not only the predicted ozone depletion, but also the accu- 
racy of this prediction. 

There are two kinds of possible inaccuracy in the predictions of ozone deple- 
tion from computer models. One kind may arise if the model fails to include 
some physically significant process, for instance, an important chemical re- 
action. Possible inaccuracies of this kind cannot be evaluated a priori, but 
can be identified by careful comparison of model predictions with actual atmo- 
spheric measurements. 

A second kind of potential inaccuracy, however, is amenable to a priori eval- 
uation. This kind of inaccuracy arises because the model input parameters are 
experimentally determined quantities, each of which has an uncertainty asso- 
ciated with it. The propagation of tliese various experimental uncertainties 
through the model calculation leads to a net uncertainty In the model outputs. 

In prhiciple, if tlie distribution of uncertainty can be estimated for each of the 
input parameters, then standard uncertainty propagation techniques could be 
used to estimate the uncertainty distribution of the output. In practice, however, 
the propagation of uncertainties through a model calculation is difficult because 
of the non-linear feedbacks involved in the radiation absorption in the strato- 
sphere. The magnitude of some of the uncertainties precludes the use of stan- 
dard ex’ror propagation formulas which are derived assuming small errors. A 
direct way to surmount these difficulties is to use a Monte-Carlo analysis, in 
which the input variables can be randomly chosen according to Uieir individual 
probabiliiy distributions. The model must then be rim enough times to generate 
a statistically valid distribution in the calculated output. This approach neces- 
sitates the development of a model which is fast enough to be run many times 
without using excess computer resources. At the same time, the model must 
be sufficiently complete to simulate the complex stratospheric processes not 
only when the input parameters are at their nominal values, but also when the 
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input parameters are at their nominal values, but also when they are one or two 
standard deviations away. 

To meet these criteria, a steady-state one-dimensional stratospheric model has 
been developed using iteration to achieve self-consistency amongst the concen- 
trations of the various species. Liberal use has been made of photochemical 
equilibrium assumptions. A rather complete chemistry is used for the odd oxy- 
gen odd nitrogen (NO^), odd hydrogen and odd chlorine (Clj^) systems 

and the radiative feedback processes which work through the dissociation coef- 
ficients are included. This paper (paper I) describes the model, while paper II 
(Stolarski, et. al. , 1977) discusses the use of this model for uncertainty analysis. 

DESCRIPTION OF MODEL 


The model developed for uncertainty studies covers the altitude range from 15 
to 50 km in 5 km steps. Because of the use of iteration to achieve self-consis- 
tency only single uncoupled diffusion equations are solved. The solutions for 5 
km altitude steps have been tested and found to agree with diffusion solutions for 
2 km and 1 km steps. The basic chemistry of the odd oxygen, nitrogen, hydro- 
gen, and chlorine compounds is included together with the methane oxidation 
chain. Fixed profiles are used for the stable molecules Nj , O, , H^O, and 
and for temperatures (see Table 1 ). N,0 and CH^ concentrations are determined 
from a self-consistent diffusion calculation with fixed lower boundary conditions 
of 300 ppbv NoO and 1.5 ppmv CII 4 . Except for the species NO, NO 2 , CIO and 
CIONOt, day-night effects are suppressed by diurnal averaging. The method by 
which these exceptions are treated is given below. Reactions used in the model 
and their rates, together with estimated uncertainties, ax’e siunmarized in Table 
2 . 


Photodissociation Coefficients 


The calculation of photodissociation coefficients uses the solar flux values re- 
ported by Ackerman (1971) as modified by the measurements of Simon (1974). 

The spectrum is divided into 25 wavelength intervals between 1800 and 3200A. 
Dissociation coefficients for the wavelength region be 3 mnd 3200A are calculated 
assuming no attenuation of the incident solar radiation. The absorption of the 
incident radiation below 3200A is assumed to be controlled by O, and O 3 while 
all other species are assumed optically thin. Table 3 shows the cross sections 
used as a function of wavelength; these values can be readily used along with spe- 
cies concentrations given later in this paper to demonstrate the accuracy of this 
assumption. 

The first two wavelength intervals from 1800 to 2000A contain the Schumann- 


Runge absorption bands of O,. These give rise to a temperafure dependent 
absorption coefficient as described by Hudson and Mahle (1972). Their technique 
was used to calculate detailed dissociation rates throughout this region for dif- 
fering Oj and O 3 column contents as a function of temperature. These were then 
combined into lOOA intervals and best-fit effective cross-sections were obtained 
as a function of temperature. Because the model extends to only 50 km, this is 
sufficient detail for all significant processes. A reexamination of these approxi- 
mations would be necessary to extend the model much beyond 60 km. 

The calculated dissociation coefficients were diurnally averaged for the spring- 
fall season at 30° latitude by the technique described by Rundel (1977). Briefly 
this consists of recognizing that in a given wavelength interval the ratio, R, of 
the attenuated solar intensity to the solar intensity at the top of the atmosphere 
is a function only of the optical depth. This function of optical depth can be di- 
urnally averaged numerically yielding an empirical formula for use in the model. 
At each altitude and within each wavelength band the optical deptlis of 0, and 
O 3 are calculated and the appropinate diurnally averaged value of R is computed. 
The dissociation (J) J coefficients for each species and each wavelength are then 
calculated and the integration over wavelength carried out. This results in con- 
siderable time saving over the technique of calculating the wavelength integrated 
J coefficient and then averaging that quantity over time of day. Figure 1 shows 
the calculated J coefficients as a function of altitude together with the standard 
case ozone profile. These values represent 24 hour averages without any effects 
due to Rayleigh scattering or sux'face x'eflection (see e.g. Luther and Gelinas, 
1976). 

Odd Oxygen 


The term "odd oxygen" is customarily used to refer to the sum of the O and O 3 
concentrations. Odd oxygen is produced by the photolysis of Oj ; it is destroyed 
by the reaction of 0 with O 3 and by catalytic cycles involving Cl^., NO^^, and HO^^. 
The calculation of destruction rates for the 0+0 3 reaction and the NOj^ and Cty 
cycles is straightforward, but for HO^^ it is somewhat complex. For tills reason, 
a generalized concept of odd oxygen has been used and is discussed below in the 
HOj^ section. This concept simplifies the determination of odd oxygen desti'uc- 
tion rates. Odd oxygen is essentially in photo-chemical equilibrium above 28 km. 
Below this altitude, chemical time constants are long compared to transport 
times and the equilibrium appro.ximation is inappropriate. Flowever, the produc- 
tion and loss rates of odd oxj'^gen have altitude distributions in this region which 
coincidentally give reasonable ozone distributions. Invoking the photochemical 
equilibrium assumption results in a large saving in computer time which in turn 
makes die Monte-Carlo uncertainty analysis feasible. 
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Figure 2 shows the production and loss rates of odd oxygen due to various reac- 
tions. It is interesting to note the differences between this breakdown of odd ox- 
ygen loss processes and those calculated 2 or 3 years ago by most modelers. 
Representative of circa 1974 model results are those of McElroy et al. (1974) 
which show the odd oxygen loss rate dominated by NOj, reactions between 30 and 
40 km by about a factor of 5. Around 45-50 km 0+0^ dominates and reac- 
tions don't become dominant until about 55 fcn. The present results indicate a 
niich larger role for mainly because of the I’educed rate now used for the 
reaction of OH and HO,, a rate which is supported by the OH measurements of 
Anderson (1976). This larger role for HO^^ reactions leads to lower calculated 
ozone concentrations around 50 km. These concentrations are in conflict with 
existing data such as that summarized by Krueger and Minzner (1976) and more 
particularly with the recent stellar occultation measurements of Riegler et. al . 
(1977). This problem has been recognized by a number of modeling groups (see 
e.g. Frederick (1977), Liu (1977), Sze (1977), Wofsy (1977)). 

Odd Nitrogen 


The concentration of odd nitrogen (N, NO, NO,, NOj, CIONO,, HNO^) is crucial 
because the NO-NO, catalytic cycle (reactions 14 and 19 in Table 2) is the major 
ozone destruction mechanism in tlie natural stratosphei-e. Odd nitrogen is pro- 
duced in die stratosphere by the reaction of 0('D) with N 2 O which diffuses up- 
ward from the earth's surface and in the mesosphere by the effects of energetic 
particles or ionizing radiation with subsequent ion- molecule reactions. It is de- 
stroyed by diffusion to the troposphere followed by rainout, and by photolysis of 
NO followed by the reaction of N('*S) atoms with NO to reform N 2 in the 35 to 65 
km region. The available data, shown in figure 3, on total odd nitrogen is sparse 
but increasing rapidly. Above 45 km, measurements of NO are equivalent to to- 
tal odd, nitrogen, since existing photochemical models, including the present 
one, predict NO to odd nitrogen ratios greater than . 9 in the daytime. Between 
35 and 45 km measurements of NO and NO, are sufficient to define odd nitrogen. 
At lower altitudes simultaneous measurements of NO, NO, , and HNO^ are re- 
quired to define the total odd nitrogen coiicentration. Additionally, at higher lat- 
itudes the amount of NtOj may be significant although it has not yet been mea- 
sured and is not included in this model. As can be seen from figure 3, there is 
a significant downward fliix of odd niti’ogen above 60 km from the mesopheric 
source. The data indicates that by 50 km this flux has reversed sign and odd ni- 
trogen is flowing upward from the stratospheric source to the sink due to photol- 
ysis of NO. The upper boundary in the present model is at 50 km, and the bound- 
condition assumed there is an upward flux of 5 x 10^ cm‘^ sec'* which has been 
adjusted to obtain a reasonable fit to the data. 

Below 50 km, the loss rate of odd nitrogen is calculated by assuming that N(‘*S) 
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atoms created by photoli of NO are m photochemical equilibrium. The loss 
processes for N(‘*S) are reaction with NO (reaction 33) converting two odd nitro- 
gens back to N 2 and reaction with O 2 or O 3 (reactions 34 and 35) reforming NO. 
Since N(‘’S) and NO exist only in the daytime, it is necessary to use daytime av- 
erages instead of diurnal averages in calculating the loss rate oi odd nitrogen. 
This has been accomplished in die present model by using the diurnal average 
results and noting that the NO/NO, ratio in the daytime is twice the diurnal av- 
erage. 

Calculation of the dissociation coefficient for NO is extremely difficult because 
absorption occurs in predissociated bands which overlap the Shumann-Runge 
bands of O 2 . The dissociation rate is thus dependent on the details of the overlap 
of these two band systems. The most complete study to date is that of Cieslik 
and Nicolet (1973) and their results at 60° solar zenith angle have been used as 
representative of the daytmie average. In this case, the diurnal averaging tech- 
nique of Rundel (1977) is not used because it does not apply to non-Beer’s law 
absorption. 

The lower boundary condition or odd nitrogen is a fixed mixing ratio at 15 Ion of 
1 ppbv (see Evans, et al. , 1976). 

Odd Hydrogen 


Odd hydrogen (H, OH, HC^ , 2 x H 2 ^ 2 ’ HNO^) is produced by reaction of 

0*(D) with H 2 O (reaction 27), Hj (reaction 28), and CH 4 (reactions 29 and 30), 
by reaction of Cl with (reaction 4) and CII 4 (reaction 6 ) and by either photol- 
ysis of or reaction of O with CH,0 to produce HCO. Although its chemistry is 
complex, the treatment of odd hydrogen is not computationally difficult in a mod- 
el because it reacts with itself to reform H 2 O rapidly enough to justify the as- 
sumption of photochemical equilibrium. The major lo jS reactions are OH + IIO 2 
and OH + ^^ 2^2 (i^eactions 39 and 40). The chemical lifetime of HjO^ in the lower 
stratosphere (15-20 km) is comparable to the transport time and the photochem- 
ical equilibrium assumption is only approximate at that level. In the present 
model CH is selfconsistently diffused, and the chemistry cf its oxidation pro- 
ducts is treated in detail. Fixed profiles of H^O and H^ assumed. 

The major calculational problem involving the odd hydrogen family is the deter- 
mination of the odd oxygen destruction rate. For odd nitrogen, this is relatively 
easy because there is only one significant cycle and its rate determining step is 
readily apparent. For odd hydrogen, however, tliere is an intricate web of inter- 
woven catalytic destruction cycles. One possible solution is to avoid the use of 
the concept of a catalytic cycle and count one odd oxygen destroyed for each of 
the reactions. 
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O + HO, ^ OH + O, 

O, + HO^ ^ OH + 20 ^ 

OH + 0 ^ O2 + H 
OH + O3 -> HO2 + O2 
H + 03-^ OH + O2 

However this neglects a reaction sequence pointed out by Crutzen (1973) which 
is important in the lower stratosphere and troposphere, namely 

NO + HO, NOj + OH 

followed by 

NO, + Iw -* NO + 0 

which produces odd oxygen and reduces the catalytic effect. A simple scheme is 
presented here to circumvent all such difficulties. The scheme consists of gen- 
eralizing the concept of odd oxygen to include more species than just O and O 3 . 
When an O, molecule is dissociated the O atoms formed may follow many paths. 
Tor instance, one may attach to O, to form O 3 . The O 3 at some point may react 
with NO, but in this reaction the oxygen atom is simply transferred to NO form- 
ing NO 2 . Photolysis of NO^ then returns the oxygen atom to its original form. 

If, however, reaction occurs between the NO^ and another oxygen atom the O 
atom on the NO^ is recombined with the other O atom for a net loss of two odd 
oxygens. The identification of the rate determining step of this catalytic cycle 
naturally follows from considering the extra O atom of NO^ as an odd oxj'gen. 
This scheme is readily generalized by assigning an odd oxygen number to each 
molecule and then simply balancing each chemical equation to determine if odd 
oxygen is produced or lost in the reaction. Table 4 shows one reasonable set of 
odd oxygen number assignments. Note that primarily short-lived radicals are 
considered odd oxygen and tha<- the total concentration throughout the strato- 
sphere. Column 4 of Table 1 indicates the change in odd oxygen number for each 
reaction used in the model. 

Odd Chlorine 


The principle reason for the development of this model is the evaluation of un- 
certainties in model predictions of chlorine perturbations. To accomplish this, 
the model is iterated to steady state with an ambient chlorine mixing ratio pro- 
file; then, additional chlorine sources are included and the model is iterated to 
a new steady state. Source calculations for Cl^^ are performed by solving the 
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( 'ffusion equation for each source molecule assuming a fixed mixing ratio at ' 
km and zero flux at 50 km. The source molecules originate in the troposphere 
and are destroyed in the stratosphere by photolysis and for CH^Cl by reaction 
with OH. They are assumed to decompose completely at the altitude where they 
initially react. The Cl^ distribution is calculated by solving the diffusion equa- 
tion with this source, a lower boundary condition of .1 ppbv and zero flux as an 
upper boundary condition. For the ambient case, CCl^ (.1 ppbv) and CH^Cl (1.0 
ppbv) are assumed to be the source molecules yielding an asymptotic Clj^ mixing 
ratio of 1.5 ppbv. The perturbing source molecules are F-ll and F-12 in the 
relative amounts indicated by the 1975 release rate data (45% F-ll and 55% F-12; 
NRC Panel on Atmospheric Chemistry, 1976) and in total amounts varying to give 
different Cl^^ perturbations. 

Chemical reactions of Cl, CIO, and HCl, as well as ClONO^suggested by Row- 
land et al. (1976), are considered. Day-night effects must be included in detail 
for an accurate calculation of chlorine nitrate concentration. This is accom- 
plished by (1) assuming that the concentration of HCl does not vary on a 24 hour 
time scale, (2) approximating the day as twelve hours of constant insolation and 
twelve hours of darkness, and (3) assuming rapid equilibration of the Cl/ClO 
ratio. These assumptions permit analytic expressions to be written for the time- 
dependent day and night CIONO, concentrations. These expressions are assumed 
equal at sunrise and sunset and sepai’ately averaged over twelve hours. Corre- 
sponding CIO day and night average concentrations are obtained, and the daytime 
value is used for reactions in which the other reactant specie is present only in 
the daytime (e.g. O or NO). OCIO is not considered because its rate of forma- 
tion by ClO+Oj seems to be too slow to make it an important species (Rundel and 
Stolarski, 1976). ClOO formation will not significantly effect model results and 
is not calculated. For a discussion of the unplications of the chemistry of OCIO 
and ClOO see Stolarski and Cicerone (1976). 

MODEL RESULTS 


Ambient Atmosphere 


Using the reaction set of Table 2 and the boundary conditions and approximations 
described above, the model was iterated until steady-state was achieved. Figure 
4 shows the resulting m>xing ratio profiles for CH^ , N^O, and CO, and the as- 
sumed profiles for Hj and H 2 O, The eddy coefficient used was two times that of 
Hunten (see Johnston, et al. , 1976), rising cxijonentially from 4 x 10^ at 15 km 
to 2 x 10^ at 50 km. The resulting mixing ratio of CH^ at 50 1cm is 200 ppbv, in 
reasonable agreement with the measurements of Ehhalt, et al. (1972, 1974). 

The calculated N 2 O profile is in die range of the mid- latitude results of Schmel- 
tekpof, et. al. , (1977). 


7 


Figure 5 shows the calculated number densities of tlie HO^ constituents H, OH, 
HOj, and H,Oj plus those of O and 0(^D). Figure 6 shows the densities of odd 
nitrogen constituents, NO, NO, and HNO 3 together with the total odd nitrogen 
mixing ratio. Figure 7 shows the I-ICl, CIO and Cl when the Cl^^ asymptotic mix- 
ing ratio is 1 . 5 ppbv. The odd chlorine mixing ratio profile is also shown. 

Effects of Chlorine Perturbations 

Calculations of the stratospheric Cl^ profile indicate onlj* minor differences in 
the profile shape for different sources (see e.g. Cicerone et al. , 1975), This 
profile begins just above die tropopause with a very low value of .1 ppbv, in- 
creases rapidly and then levels off to a constant asymptotic mixing ratio. Mea- 
surements of HCl by Lazurus, et al, (1975), Williams et al. (1976), Farmer^ 
al. , (1976), Ackerman et al. (1976) and Eyre and Roscoe (1977) tend to confirm 
this shape. One of the profiles of Lazrus (1977) shows a falloff at high altitudes 
larger than can be explamedby conversion to CIO, Additionally the concentration 
of CIO measured by Anderson et al. (1977) is larger than would be expected ei- 
ther from known sources or from the measured HCl concentration. These mea- 
surements arc thus inconclusive so until more frequent and simultaneous obser- 
vations of all chlorine species are made, it seems reasonable to assume the 
shape described above to be correct. Accordingly, runs of the model were made 
with an ambient Cl^^ profile, and then with varying amounts of chlorine added. 

The mechanism for this was always an F -11 and F-12 addition at the lower bound- 
ary of the source calculation in the ratio 45% to 55%. The change in the ozone 
content due to the addition of chlorine, is shown in Figure 8 as a function of the 
asymptotic CIX mising ratio for a standard case and for models where the as- 
sumed H^O profile or the N^O lower bcundar 3 ' condition was varied by a factor 
of 2. Although these .suits were generated for F-li and F-12 they should be 
relatively independent of the source. The s >pe of the curve in figure 8 can be 
interpreted as the ozone destruction efficiency. This concept has been previously 
used effectively by Wofsy and McElroy (1974) and Donahue et al. (1976). The re- 
sults for large chlorine perturbations are shown as a model diagnostic. Any re- 
sults with more than about 10 % ozone depletion should be viewed with skepticism 
because of the possibility of a breakdown of model assumptions. 

Effects of NO^ Perturbations 

The model was also run with an additional injection of NOj^ at 20 km to simulate 
the potential effect of an SST fleet. Because of the altitude step size, the injec- 
tion is in a 5 km wide band around 20 km. Figure 9 shows the results for the 
same five cases. The curve crossings at low NO^ injection are caused by CIONO 2 
formation decreasing the catalytic efficiency of the ambient Cl,,. 

A. 
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SUMMARY AND CONCLUSION 


This paper has described a one-dimensional, steady-state, diurnally-averaged 
model of stratospheric photochemistry specifically designed to address the ques- 
tion of the propagation pf reaction rate uncertainties and the implications on the 
resultant precision of model results. This model has been shown, despite its 
simplicity and calculational speed, to adequately reproduce the effects calcu- 
lated using significantly more complex models. Thus an uncertainty propagation 
analysis performed on this model should simulate reasonably well the results 
which would be obtained by such an analysis on a more complex model. The 
analysis will probably never be carried out on a "complete” model because of 
computer time limitations. 

Acknowledgements 


We gratefully acknowledge useful discussions with Robert D. Hudson, Robert T. 
Watson and Ralph J. Cicerone. We thank Donald J. Kessler who determined the 
temperature dependent curve fits for tlie Schumann-Runge Bands. One of us 
(DMB) wishes to thank the National Academy of Sciences for a Resident Research 
Associateship at Goddard Space Flight Center during which most of this work was 
carried out. 


9 


REFERENCES 


1. Aclverman, M. , "Ultraviolet solar radiation related to mesospheric pro- 
cesses" in Mesospheric Models and Related Experiments ed. by G. Fiocco, 
D. Reidel, piibl, Dordrecht, Holland p. 149- (1971). 

2. Ackerman, M. , J. -C. , Fontanella D. Frimout, and A. Girard, "Simul- 
taneous measurements of NO and NO^ in the stratosphere", Planet. Space 
Sci., 651-660, (1975). 

3. Ackerman, M. , D. Frimout, A. Girard, M. Gottignies and C. Muller 

"Stratospheric NCI from Infrared Spectra", Geophys. Res. Lett. , 81 

(1976). 

4. Anderson, J. G. , J. J. Margitorn and D. H. Stedman, "Free Chlorine in 
the Stratosphere: An IN SITU Study of Cl and CIO", unpublished manu- 
script, Febi’uary, 1977. 

5. Cicerone, R. J. , D. H. Stedman and R. S. Stolarski, "Estimate of late 
1974 stratospheric concentration of gaseous chlorine compounds (CIX)", 
Geophys. Res. Lett. , 2, 219 (1975). 

6. Crutzen, P. , "A discussion of the chemistry of some minor constituents 
in the stratosphere and troposphere", PAGEOPH 106-108, 1385-1399 
(1973). 

7. Donahue, T. M. , R. J. Cicerone, S. C. Liu, and W. L. Chameides, 
"Effect of add hydrogen on ozone depletion by chlorine reactions", Geo- 
phys. Res. Lett., £, 105-108, 1976. 

8. Drummond, J. W. , J. M. Rosen, D. J. Hofmann, "Balloon-borne chem- 
iluminescent measurement of NO to 45 km", Nature, 265 , 319 (1977). 

9. Ehhalt, D. H. , L. E. Heidt, and E. A. Martell, "The concentration of 
methane between 44 and 62 kilometers altitude", J. Geophys. Res. , 77, 
2193- (1972). 

10. EMialt, D. H. , L. E. Heidt, R. H. Lueb, and N. Roper, "Vertical pro- 
files of CH^, Hj, CO N,0 and CO, in the stratosphere", Proc. Third 
Confei'ence on CIAP , ed. by A. J. Broderick and T. M. Hard. DOT-TSC- 
OST-74-15, pg. 153-160 (1974). 

11. Evans, W. F. , J. B. Kerr, and D. I. Wardle, "The AES stratospheric 


ballon measurements project: preliminary results", Proc. Fourth Confer- 
ence on CIAP, ed. , by T. M, Hard and A, J. Broderic, DOT-TSC-OST- 
75-38, pg. 412-416 (1976). 

12. Eyre, J. R. and H. K. Roscoe, "Radiometric measurement of strato- 
sphei’ic HO", Nature, 266 , 243-244 (1977). 


13. Framer, C. B. , O. F, Raper and R. H. Norton, "Stectroscopic detection 
and vertical distribution of HCl in the troposphere and stratosphere", Geo- 
phys. Res. Lett. , 3, 13 (1976). 

14. Frederick, J. E. , P. B. Hays and S. C. Liu, "The need for a revised 
theory of mesospheric ozine", EOS, in press (1977). 

15. Horvath, J. J. and C. J. Mason, personal communication (1976). 

16. Hudson, R. D. and S. H. Mahle, "Photodissociation rates of molecular 
oxygen in the mesosphere and lower thermosphere", J. Geophys. Res. , 

77, 2902 (1972), 

17. Johnston, H. S. , D. Kattenhorn and G. Whitten, "Use of excess carbon 
14 data to calibrate models of stratospheric ozone depletion by supersonic 
transports", J. Geophys. Res., 368-380 (1976), 

18. Krueger, A. J. , and R. A. Minzner, "A mid-latitude ozone model for the 
1976 U.S. standard atmosphere", J. Geophys. Res. , W, 4477 (1976). 

19. Liu, S. C. and R. J. Cicerone, "Comparison of theory and measurement 
of ozone above 35 km", EOS in press, 1977. 

20. Lazrus, AIL., B. W. Gandrud, R. N. Woodard, and W. A. Sedlacek, 
"Stratospheric Halogen Mea.surements", Geophys. Res. Lett. , ^39 
(1975). 

21. Lazrus, A. , Position paper on Hydrogen Halide Compounds in the Strato- 
sphere presented at NASA CFM assessment workshop, Warrenton, Va. , 
Jan. (1977). 

22. Luther, F, M. and R. J. Gelinas, "Effect of molecular multiple scatter- 
ing and surface albedo on atmos"'heric photOdissociation rates", J. Geo- 
phys. Res., 1125-1132 (1976). 

23. McElroy, M. B. , S. C. Wofsy, J. E. Penner and J. C. McConnell, 


11 


"Atmospheric Ozone: Possible Impact of Stratospheric Aviation", J. Atmos 
Sci. , 81, 287-303 (1974), 

24. Meira, L. , "Rocket measurements of upper atmosphere nitric oxide and 
their consequences to the lower ionosphere", J. Geophys. Res., 202 
(1971). 

25. National Research Council Panel on Atmospheric Chemistry. , Halocarbons: 
Effects on Stratospheric Ozone, Washington, D.C., September, 1976. 

26. Riegler, G. R. , S. K. Atreya, T. AI. Donaliue, S. C. Liu, B. Wasser and 
J. F. Orake, "UV Stellar Occultation Measurements of Nighttime Equa- 
torial Ozone", Geophys. Res. Lett., in press (1977). 

27. Rowland, F. S. , J. E. Spencer, and M. J, Molina, "Stratospheric Forma- 
tion and Photolysis of Chlorine Nitrate", J. Phys. Ciiem. , 2711-2713 

(1976). 

28. Rundel, R. D. and R, S. Stolarski, "A reexamination of the photochemistry 

of the Cl, - O^ system", J. Geophys. Res. , 5759-5764 (1976). 

29. Rundel, R. D. , "Determination of diurnal average photodissociation rates", 
J. Atmos. Sci. , (1977), in press. 

30. Schmeltekopf, A. L. , D. L. Albritton, P. J. Crutzen, P. D, Goldan, W. 

J. Harrop, W. R. Henderson, J. R. McAfee, M. McFarland, II. I. 

Schiff, T. L. Thompson, D. J. Hofmann and N. T. Kjome, 'Btratospheric 
Nitrous Oxide Altitude Profiles at Various Latitudes", Position paper pre- 
sented at NASA CFM assessment workshop, Warrenton, VA, January 
(1977). 

31. Simon, P. , "Balloon measurements of solar fluxes between 1960 A and 
2300 A", in Proc. Tltird Conference on CIAP , ed, by A. J, Broderick and 
T. M. Hard, DOT-TSC-OST- 74-15, O. 137-142 (1974). 

32. Stolarski, R. S. and R. J, Cicerone, "Chlorine in the stratosphere", Proc. 
Fourth Conference on CIAP , ed. by T. M. Hard and A. J. Broderick 
DOT-TSC-OST- 75-38, pg. 280-285 (1976). 

33. Stolarski, R. S., D. M. Butler and R. D. Rundel, "Monte-Carlo analysis 
of uncertainty propagation in a stratospheric model: II. Application to reac- 
tion rate uncertainties", submitted to J. Geophys. Res. (1977). 


12 




34. Szo, N. D»| 't:\tmospliereozone; Comparison of theory and observations*’, 
EOS, in jiveas (1977). 

35. Tlsone, G. , ^■'Measurements of NO densities during sunrise at Kauai*', J. 
Geophys. NeS. , *" y 746- (1973). 

36. Wofsyii Sk C. a^ud M, B. McEiroy, "HOj^, NOj^and ClOj^: their role in 

atntospherio photoshenilstry", Can. J. Chem. , 1582 (1974), 

37i Wofsy, S, C. , "Temporal and latitudinal variations of stratospheric trace 
gases: A critical comparison between theory and experiment*', submitted 
to J. Geophys. Res. (1977). 


13 



Table 1 


FiAed Input Profiles 


Altitude 

(kni) 

N2 

(cm'^) 

°2 

(cm'^) 

H,0 

(cnT^) 

(cnT^) 

T 

(°K) 

15 

3.20(18) 

8.49(17) 

1.28(13) 

1.99(12) 

215 

20 

1.46(18) 

3.87(17) 

6.81(12) 

8.70(11) 

216 

25 

6.60(17) 

1.75(17) 

3.53(12) 

3.69(11) 

221 

30 

3.03(17) 

8.02(16) 

1.84(12) 

1.59(11) 

226 

35 

1.39(17) 

3.69(16) 

9.31(11) 

7.12(10) 

236 

40 

3.24(16) 

8.59(15) 

2.34(11) 

2.32(10) 

264 

50 

1.69(16) 

4.48(15) 

1.18(1 1) 

1.78(10) 

271 
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Table 2 

Reaction Rates and Uncertainties 



REACTION 

RATE 

(cm^ sec’* or cm®sec’*) 

UNCERTAINTY 
(la in log k) 

A Odd 
Oxygen 

1. 

a+ 03 ->ciio +02 

2.7-10-**e-2^’/'^ 

.06 

0 

2. 

CIO+O-^Cl+Oj 

7.7-10’** 

.06 

-2 

3. 

aiO+NCHCl+NOj 

2.2-10’^* 

.15 

0 

4. 

Cl+H,-^HCi+H 

3.5-10’“ 

.04 

0 

5. 

CI+HO^-^HCl+Oj 

3.0-10’** 

.30 

0 

6. 

C1+CH^->HCI+CH3 

7.3.10-12 e- 1260 /'^ 

(+).06 

(-).18 

0 

7. 

CL+H202->'HC1+H02 

1.7-10-12e-384/T 

.18 

0 

8. 

HCl+OlHCl+H^O 

3 . 0 . 10 - 12 c-'* 25 /T 

.05 

■1 

9 . 

ClO+NOj+M-^ClONOj+M 

5 .MO’”e^l 030 /T 

(+).06 

(-).15 

0 

10. 

ClONO^+O^CIO+NOj 

4.5.10-' 2 

.30 

0 

11. 

ClONOj+Oll-^products 

1.2.10-12 e- 333 /T 

.18 

0 

12. 

O+O^+M-^Oj+M 

1.07.10-2'’ 

.07 

0 

13, 

O+O3-O2+O2 

I. 9 .IO-" C- 2200 /T 

.10 

-2 

14. 

NO+Oj-i-NOj+Oj 

2.M0-'2e-'‘’20/1’ 

.08 

0 

15. 

N0,+0,->-N0,+0, 

1.2.10-12 e 2450 /T 

.04 

0 

16. 

OH+O3-+HO2+O2 

1.5.10-12 e-llll/T 

.15 

-2 

17 . 

H0,+0,-*0H+0,+0, 

1.0.10’‘2 e-'525/T 

.30 

0 

I''. 

H+O^-^OH+Oj 

1.2.10-11 e-2110/T 

.15 

0 

19. 

NO., +0->'N0+02 

9.1.10-12 

.03 

-2 

20. 

N 02 + 0 +M->N 03 +M 

1.0.10"21 

.20 

0 

21. 

OH+O-^-H+Oj 

l.O-lO-lle-220/T 

.15 

-2 

22. 

HO^+O-^OH+Oj 

1.0.10-11 e-220/T 

.15 

0 

23. 

H^Oj+O-KDH+HOj 

2.75.10-12 ^- 2125 /T 

.15 

0 

24. 

HN 03 + 0 ->OH+N 03 

i.o.io'i'i 

(+).15 
(.).!. 0 

0 

25. 

CH^O+CKHCO+OH 

2.0.10-11 e-l'’20/2- 

.20 

0 

26. 

N20+0(^D)^NO+NO 

5 . 5 . 10-11 

.10 

+ 1 
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Table 2 Reaction Rates and Uncertainties (Continued) 



REACTION 

RATE 

(cm^scc”' or cni^ sec'^) 

UNCERTAINTY 
(Iff in log k) 

— 

27. 

H20+0 (‘d)^0 H+0H 

2.3- 10'^” 

.05 

■H 

28. 

H2+0 (‘d)-^OH+H 

9.9-10"^‘ 

.05 

0 

29. 

CH^+0(’ D)->0H+CH3 

1.3-10“^'’ 

.05 

0 

30. 

CH^+ 0 (‘D)^H 2 +C 1 - 1,0 

1.4-10'^’ 

.10 

0 

31. 

0(^D)+M->0(^P)+M 

2 . 0 - 10 -' 

.05 

0 

32. 

NO+NOj^NO^+NO., 

8.7-10"'^ 

.70 

0 

33. 

N+NO^Nj+O 

8 . 2 - 10 -" c-'""/'!' 

.10 

0 

34. 

N+0.,->-N0+0 

5 5 . 10-12 ^-3220/T 

.10 

-t -2 

35. 

N+0j->N0+0, 

5.0-10-'2e-®^“/''' 

.30 

0 

36. 

NO+HO^-^NO^+OH 

1.5-10-“ 

.20 

+2 

37. 

NO^+OH+M-^HNOj+M 

(see expression below) 

.05 

0 

38. 

OH+HNO 3 -+NO 3 +H 2 O 

8 . 0 - 1 0 "'*' 

.05 

-1 

39. 

OH+H^O^^HO^+H^O 

1.0-10-“ 

.15 

-1 

40. 

0 H+H 02 -> 02 +H 20 

3.0-10"“ 

.25 

-1 

41. 

OH+OIHO+HjO 

1 . 0 - 10 -“ e-^^°/’'' 

.20 

-1 

42. 

oh+M2->h+h.,o 

8 . 0 - 10 -'^- 2 ' 0 °/^ 

.10 

-1 

43. 

OH+CH^-^CHj+H^O 

2.36-10-' 2 

.10 

-1 

44. 

OH+CH^O^HCO+H^O 

3.0-10-' ' 

.20 

-1 

45. 

0 H+Cll 3 CKCH,Cl+H ,0 

2_2.io-*2e-‘142/T 

.10 

-1 

46. 

OH+CCKCOj+II 

1.4-10-'^ 

.10 

-2 

47. 

OH+OH+M-^-HjO^+M 

1.25-10-'^^^°°/'^ 

.15 

-2 

48. 

H+O^+M-^HOj+M 

2 . 1 - 10 -'- e^-'"/''' 

.10 

0 

49. 

I 10 ,+l-I 0 .,->ll., 02+02 

5.0-10-"e-^‘'°/'’ 

.30 

0 

50. 

CW 3 + 02 +M->CH 30 ,+M 

3.85-10"" 

.50 

+2 

51. 

CI-l 30 ,+N 0 ^CU 30 +N 02 

3.3-lQ-'^ 0"'°°^''' 

1.0 

0 

52. 

CI-l 302 +N 02 ->CH 30 +N 03 

1.0-10'" 

1.0 

0 

53. 

CH 302 +H 02 ->Cli 300 H +0 

2 6.7-10"''' 

1.0 

0 

54. 

CH 30 + 02 ->CH 20 +HO, 

16.10-13, -3300/T 

1.0 

0 

55. 

HCO+02^CO+H02 

6.0-10"'- 

.20 

0 
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Table 2 Reaction Rates and Uncertainties (Continued) 


for HNOj formation the following expression was employed: 
logjQk = -AT/(B+T) 

where: A = -10.4432 + 6.33Z- ,4321 26Z^ +8.47677 -10‘^Z^ 
B = 327.372 + 44.5586Z - 1.38092Z" 

Z = logjg [M] 
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Tabic 3: Solar Flux and Crossections 

Solar (lux is in photons/cm^ sec; Crossection arc in cm^; Visible J Coeffident are lii sec '. 






















































































































Table 4 



Odd 

Odd 

Odd 

Odd 


Oxygen 

Nitrogen 

Hydrogen 

CliJorine 


1 

0 

0 

0 

0 

1 

0 

0 

0 

N 

0 

1 

0 

0 

NO 

1 

! 

0 

0 

NOj 

2 

1 

0 

0 

NO3 

3 

1 

0 

0 

HNO3 

3 

1 

I 

0 

H 

0 

0 

1 

0 

OH 

1 

0 

1 

0 

HO3 

0 

0 

1 

0 

»2°2 

1 

0 

0 

0 

CO 

1 

0 

0 

0 

CH3 

0 

0 

1 

0 

CH303 

2 

0 

I 

0 

CH30 

1 

0 

I 

0 

H2CO 

1 

0 

0 

0 

HCO 

1 

0 

1 

0 

CH3OOH 

2 

0 

2 

0 

Cl 

0 

0 

0 

1 

HCl 

0 

0 

1 

1 

CIO 

I 

0 

0 

1 

N2O3 

5 

2 

0 

0 

ClONOj 

3 

I 

0 

] 


I The following species jii the model have only zero odd species numbers: 
■ Oj, Nj, COj, CH^, NjO, HjO, 


f ’ 

ii 
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Fi^iurc 2a. Production rates (a) and loss frequencies (b) of odd oxygen due to various processes as a function 
of altitude for Cl^ = l.fi ppbv. Loss frequency is defined as the volume loss rate in cm~^ sec~^ 
divided by the odd oxygen concentration in cm~^ . 




Figure 2b. 
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Figure 4. Mixing ratio profiles for CH^, N2O, CO, H2 and H2O. 
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DENSITY (cnr^) 

Figure 5. Calculated 24 hour average concentrations of HO^ constituents plus O and 0('D). 





Figure 6. Cal 
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Calculated 24 hour average chlorine concentrations for an asymptotic mixing ratio of 1 . 5 ppbv CIX 
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Figure 8. Fractional column ozone depletion as a function of asymptotic addc(i 
Cl)( mixing ratio. 


28 



COLUMN 



NOx PERTURBATION 
(Injected Flux cm ^ sec 
with 1.5 ppb Clx 


Figure 9. Fractional column ozone depletion as a function of added N(>x flux in 
jected in a 5 km band about 20 km. 



